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The problems of increasing the operating efficiency of vibratory mechanisms by means of multiple combina¬ 
tion parametric resonance are examined. It is proposed that this can be done with a parametric vibration ex¬ 
citer where the driving force is generated by the interaction of masses, which rock in the centrifugal inertial 
force field generated by a rotating rotor and in the field of gravity, with a working unit which elastic compo¬ 
nents secure to a base. 


Vibratory feeders and conveyors [1], vibratory mecha¬ 
nisms which initiate the flow of loose materials, and other vi¬ 
bratory mechanisms are the principal components of produc¬ 
tion-flow lines conveying the components of the glass batch 
and cullet in sectional and machine-tank factor works. 

Electrical, mechanical, and pneumatic vibratory units are 
used as vibrators in such mechanisms. Electromagnetic vi¬ 
brators, in which ac electro- or permanent magnets are used 
to produce the vibrations, and electromechanical vibrators, 
where unbalanced masses placed on the shaft of an electric 
motor are used to produce vibrations, are the vibrators most 
widely used in the glass industry. 

Most vibratory mechanisms and feeders of the glass- 
batch components operate in a forced-vibration regime at 
frequencies far from resonance, as a result of which energy is 
not used efficiently. This is because in such frequency ranges 
the vibration regime of a vibratory mechanism is insensitive 
to and essentially independent of changes in the parameters 
of the technological load (a decrease of the material level in a 
hopper or fluctuations of the moisture content, granulometric 
composition, or bulk density of the transported material). 

A novel method of parametric excitation of vibrations of 
mechanical systems (RE Patents Nos. 2072660 and 2072661), 
which makes it possible to implement a multiple combina¬ 
tion parametric resonance regime, can increase the operating 
efficiency of vibratory mechanisms. 

The properties of a combination parametric resonance 
can be analyzed for the example of a dynamic model of the 
vibratory setup shown in Fig. 1. 
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The vibratory setup (see Fig. 1 a) contains a balanced ro¬ 
tor 1 with mass m 0 , rigidly secured to a drive shaft 2, whose 
rotation axis is oriented perpendicular to the vertical plane. 
The drive shaft is set in bearings which are held in the sup¬ 
ports 3 rigidly tied to the working unit 3 with mass M 0 . Elas¬ 
tic elements 6 secure the working unit to the base 5. A 
damper 7 simulates friction from the technological load. It is 
assumed that the working unit undergoes only translational 
motion along the O y axis. The coordinate system Ax'y 'z' 
whose origin lies at the center of the rotor and whose axes 
are parallel to the corresponding axes of the stationary coor¬ 
dinate system Oxyz moves translationally relative to the lat¬ 
ter. In a position of static equilibrium the Az' axis coincides 
with the Oz axis. 

The rotor has three periodically alternating open circular 
running tracks 8 , whose centers are shifted from the rotation 
axis of the rotor by the same distances AB = l (see Fig. lb ). 
Rocking bodies (runners) 9 each with mass m are placed on 
the running tracks and are free to roll in their tracks. 




Fig. 1 . Diagram of a parametric vibratory setup. 
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This setup can be regarded as a mechanism for periodi¬ 
cally changing the inertial parameters of the vibratory system 
in time with period 2k/( o. A periodic change in the parame¬ 
ters of the system is obtained by imposing on the motion of 
the runners nonstationary (rheonomic) constraints as a result 
of the uniform rotation of the rotor with angular velocity oe>. 
The rotating rotor with the runners rocking in a centrifugal 
force field forms an inertial element (IE) of the parametric 
vibration exciter (RF Patent No. 2072661). 

The dynamics of this vibratory mechanism with the rota¬ 
tion plane of the rotor of the IE in a horizontal position has 
been studied in [2-4]. When the rotor is in a vertical posi¬ 
tion the forces due to the weights of the runners create a 
nonstationary force field. The coefficients of the force func¬ 
tion are periodic functions of time with period 2k/ co, which 
results in additional parametric excitation [5]. 

The positions of the running tracks are determined by the 
angles \\r k = co t + 2nk/N (N is the number of runners), and the 
positions of the runners are determined by the angles cp^ 
(k= 1,2,..., N ). The rocking angles cp^ of the runners and the 
displacement y of the working unit are taken as the general¬ 
ized coordinates of this mechanical system. 

Moving along the O y axis the working unit overcomes 
the elastic forces of the nonlinear springs and resistance 
forces, whose characteristics are given in the form 

F y ~cy + c l y 3 ; 

R y = by+b x y 2 y, 

where c is the total stiffness of the elastic elements, c l is the 
coefficient of nonlinearity of the elastic restoring force, b and 
b { are positive constants, and an overdot denotes differentia¬ 
tion with respect to the time t. 

The following differential equations, governing the mo¬ 
tion of the vibratory mechanism, with periodic coefficients 
were derived on the basis of Lagrange’s equations of the se¬ 
cond kind: 


<P k +2(«o +K (S ?l )q>* + v 2 co 2 sin <p t = 

-v 2 (g/l+y)c os(y k +V| I k ); 
y + 2(n + hy 2 )y+(a 2 y (\ + yy 2 )y = 

N 

Zt( co+( i ) ^ ) 2 sin (<P* )-<i>* cos (<p* )l 



where n 0 and h 0 are, respectively, the coefficients of the li¬ 
near and nonlinear damping of the runners; k = 1, 2, ..., N, 
and N= 3; v = (mp c l/J B ) 0 5 (p c is the density of the glass 
batch; J B — the reduced moment of inertia of a runner rela¬ 
tive to the axis of rolling) — the relative rocking frequency 
of the runners; g is the acceleration of gravity; n = b/2M and 
h = b { /2M are, respectively, the coefficients of linear and 
nonlinear damping of the working unit; oo y = c/M is the par¬ 
tial frequency of the intrinsic vibrations of the working unit; 
y = c { /c; p c = BC k ; M - M 0 + m 0 + Nm is the total mass of 
the entire system. 


We shall investigate the resonance properties of the sys¬ 
tem (1) for co v « co/2. 

The terms containing the second and third powers of the 
coordinates and their derivatives play the main role in the 
formation of the most important effects in the nonlinear sys¬ 
tem (1). Consequently, the trigonometric functions sin cp^ and 
cos cp^ in Eqs. (1) are replaced by two terms of their series 
expansions. 

Disturbance of the position of equilibrium and self-exci¬ 
tation of vibrations are possible in systems whose motion is 
described by differential equations with periodic coefficients. 
The existence of such phenomena and their possible applica¬ 
tions in electromagnetic systems are well known. 

For our nonlinear system it is impossible to draw a sharp 
line between forced and parametric excitations. This is be¬ 
cause a forced vibratory process in a nonlinear system can 
cause the corresponding parameters to change periodically. 
The driving force in Eqs. (1) plays the role of “oscillations of 
pumping.” On this basis we introduce the new variable u k de¬ 
fined as 


(p^ = u k + d cos \\f k , 


where d cos \j i k is a particular solution of the equation 
cp /c + v 2 co 2 cp£ = - v 2 jlxcos \\f k , corresponding to forced non¬ 
resonant oscillations of the runners; d = pv 2 (l - v 2 )~ l (D~ 2 , 
and ju = g/l(£> y 2 . 

When the working unit is in a fixed position (y = 0) the 
expression dcos \\f k is an approximate solution of the sys¬ 
tem (1), corresponding to nonresonant oscillations of the run¬ 
ners. It can be shown that in this case the center of mass of 
the system of runners describes a circle with radius cP-p c /8 in 
the direction of rotation of the rotor of the IE with respect to 
the coordinate system Ax'y'z' [4, 5]. Its rotational angular 
velocity is Mo (N= 3). As a result of the rocking of the run¬ 
ners, the rotor of the IE becomes unbalanced. The rotation of 
the center of mass of the system of runners gives rise to a 
centrifugal inertial force which is transmitted to the working 
unit of the vibratory mechanism. The harmonic components 
of the force from the rocking of all runners mutually cancel, 
except for the harmonics whose frequencies are multiples of 
Mo. It is these components that are transmitted to the work¬ 
ing unit. The harmonic with frequency Mo has the largest 
amplitude, so than the effect of the higher-order harmonics 
can be neglected. Since Aco » <$ y , to a first approximation 
the effect of the harmonic with frequency Mo on the working 
unit outside the resonance zone can be neglected and it can 
be assumed thaty = 0. Thus, the system (1) in the new vari¬ 
ables u k is in a position of equilibrium u k = 0,y = 0. 

The parametric resonances with can be appear in the sys¬ 
tem (1) arise near the frequencies 

oo = (co.+ coy)j _1 , ij = 1, 2,..., N+ 1, s= 1,2,... 

The frequencies with i=j , i.e., co = 2co^ _1 , are called 
simple resonance frequencies and the frequencies with i*j 
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are called combination resonance frequencies (co, are the re¬ 
lative partial frequencies of the characteristic vibrations (RF 
Patent No. 2072661)). 

We shall now examine the combination resonance where 
vibrations are excited at the frequencies ocq and <u 2 which are 
related with the excitation frequency co by the relation 

oo = oo 1 + co 2 , (2) 





Fig. 2. Amplitude - frequency characteristics of the working unit 
as a function of the linear damping. 


where ccq - vco, oo 9 « m , and v = 1/2, co - oo/2. 

1 ^ y y 

For such tuning the partial frequencies satisfy the condi¬ 
tions for the excitation of combination and simple parametric 
resonances simultaneously. 

We shall seek the solution of the system of equations (1) 
in the new variables, taking account of Eq. (1), in the form 

u k = a k cos (oo 1 1 + 0^), y = acos (co 2 1 + 0), k= 1, 2,..., N 

and we shall use the Bubnov - Galerkin method [6] to deter¬ 
mine unknown amplitudes a k and a, the phases 0 A , and 0, and 
the generation frequencies oo 1 and co 2 . 

Since the runners are identical, a k = a 0 ,k= 1,2,..., TV. By 
virtue of the symmetry of the IE we shall seek the solution of 
the transformed equations (1) for 0^ = 2nk/N. Using the 
Bubnov - Galerkin method we arrive at a system of alge¬ 
braic equations for a 0 , a , e, ccq, and u> 2 . The solution of this 
system was obtained numerically by iteration using Seidel’s 
method. The results of the numerical solution of the equa¬ 
tions are presented in Figs. 2 and 3 with tuning v = 1/2, 
co >; « co/2 and the following values of the parameters of the 
system: 7V=3, co y = 25 sec -! , /7 0 = 0.015, h 0 = h = 0.02, 

s = 0.03, |3 0 = 0, p=1.57 (n 0 = n 0 /(o v ; n = n/ co v ; h {) = 

h 0 /co y ; h = h /co v ; z=N ( mp c ) 2 /(2 J B M) is a coefficient that 
is proportional to the ratio of the total mass of the runners to 
the mass of the entire system; p 0 = Y l 2 ? S=co/cd v , co, = 

uq /co y , co 2 = oo 2 /co , a = a 11 . 

Studying the stability of the stationary working states of 
the mechanism, it can be shown that the dashed curves in 
Fig. 2 correspond to unstable solutions. 

Analysis of the amplitude - frequency characteristics 
showed that the resonance curves are flat and do not have the 
maximum characteristic of high-Q oscillatory systems at re¬ 
sonance of forced oscillations. Using only trigonometric 
nonlinearities, which are due to the characteristic features of 
the kinematics of the rotating IE, makes it possible to in¬ 
crease substantially the stability of the resonance regime of 
the vibrations of the working unit with linear elastic ele¬ 
ments. Increasing the coefficient n of linear damping of the 
working unit of the vibratory mechanism for small values of 
the linear damping n 0 of the runners enlarges the resonance 
zone (see Fig. 2). With a six-fold increase of the linear damp¬ 
ing the vibrational amplitude of the working unit decreases 
by only a factor of 2.5. This also attests to adequate stability 
of the resonance regime of the vibrations. The amplitude - 
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Fig. 3. Generated frequencies S] (5) and co 2 (5) versus the excita¬ 
tion frequencies: 7, 2, and 3 ) n = 0.02, 0.08, and 0.12, respectively. 



frequency characteristics of the runners are similar to those 
of the working unit. 

We note that the resonance curves in Fig. 2 were con¬ 
structed neglecting the stops limiting the rocking angles of 
the runners. It is assumed that the runners have individual 
running tracks which have a large rocking angle and are posi¬ 
tioned in parallel planes. 

The variation of the generated frequencies ocq and co 2 as a 
function of the excitation frequency is shown in Fig. 3. The 
vibrational frequency co 2 of the working unit remains almost 
constant and approximately equal to the 1 or, in dimension¬ 
less units, cd 2 « co/2. 

The following threshold condition must be satisfied for 
vibrations to exist: 

8v«w 0 
8 >--. 

1 — V 

As one can see, the threshold value of the parameter s 
depends strongly on the parameter v. For example, when the 
rocking frequency of the runners changes from v = 0.25 to 
v = 0.50 the threshold value increases by a factor of 3. Con¬ 
sequently, for v = 0.50 there arises the problem of ensuring 
self-excitation of the vibrations of the vibratory machine for 
sufficiently small s. 

In addition, dry friction in the rotor - runner system and 
in the units securing the elastic elements also results in an ex¬ 
citation threshold. However, in the present case the high-fre¬ 
quency component of the vibrations converts dry friction 
into viscous friction [7]. Outside the resonance zone two pro¬ 
cesses occur in the vibratory system of this setup: nonreso- 
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nant rocking of the runners with driving frequency cd and vi¬ 
brations of the working unit with frequency 3cd (co-vibra- 
tions). In the resonance zone the runners oscillate with fre¬ 
quency corvee (v = 0.5), i.e., with a frequency approxima¬ 
tely six times lower than the frequency of the co-vibrations 
of the working unit. Hence it is evident that under the condi¬ 
tions considered here dry friction becomes viscous friction. 
For this reason, the excitation threshold due to dry friction 
can be overcome with zero initial conditions at startup. We 
note that the vibratory activity of the runners promotes 
self-excitation of the vibratory mechanism because they are 
tuned to the simple (main) parametric and autoparametric 
resonances for any parametric excitation frequency oe>. These 
factors guarantee that the setup will always enter a resonant 
vibratory regime. 

When the mechanism is started up and reaches rotational 
frequencies in the region of instability relative to the equilib¬ 
rium position u k = 0, y = 0 (resonance zone), a multiple com¬ 
bination parametric resonance is excited. As a result of the 
rocking motion of the runners in the IE a definite phase rela¬ 
tion is established between them (0^ = 27 ik/N; k = 1, 2, 3; 
N= 3) and their common center of mass moves along a 
closed curve. It is easily shown that to a first approximation 
the center of mass of the system of runners with respect to 
the coordinate system Ax 'y r z ' with origin at the center of the 
rotor in the IE will move according to the law: 




Pc a o 


Pc a 0 
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sin co 2 1 ; 


cos cd 2 t. 
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where x' c = x c , y c — y c —y; x c and y c are the coordinates of 
the center of mass of the runners. 

It is evident that the center of mass in the primed coordi¬ 
nate system describes a circle in the direction of rotation of 
the rotor in the IE with frequency co 2 * gd/2. The radius of the 
circle along which the center of mass of the runners moves 
depends nonlinearly on their amplitude. The circular motion 
of the center of mass of the runners gives rise to an unba¬ 
lanced centrifugal inertial force which also rotates with fre¬ 
quency co in the plane of the IE. Since the working unit of the 
setup is tuned to the frequency cd v = cd 2 « cd/2, the centrifugal 
force component with frequency cd 2 will give rise to reso¬ 
nance vibrations of the working unit. It can be shown that in 
this vibrational system the vibrations are coupled with one 
another. The unbalanced centrifugal inertial force arising as a 
result of the rocking of the runners gives rise to resonance vi¬ 
brations of the working unit and the vibrations of the work¬ 
ing unit give rise to the resonance oscillations of the runners. 
This bilateral interaction gives rise to a multiple combination 
parametric resonance. The working unit undergoes almost 
subharmonic vibrations of the order of 1/2, i.e., with fre¬ 
quency cd 2 « cd/2. 


The principle of operation of the vibrational mechanism 
consists of the following. When the runners rock back and 
forth the rotor in the IE automatically becomes unbalanced, 
and therefore a centrifugal inertial force is the driving force 
of the parametric vibration exciter. However, in contrast to 
the ordinary centrifugal vibrator a parametric vibration ex¬ 
citer converts the rotational energy of the IE more efficiently 
into the mechanical vibrational energy of the working unit. 
The vibratory mechanism acquires new properties which are 
not completely obvious. The resonance zone can be ex¬ 
panded by increasing the damping by a definite amount, and 
most importantly the resonance vibrational regime for a 
high-Q vibratory system becomes much more stable. 

The fact that the frequency of the mechanical vibrations 
of the working unit of such a setup can be decreased without 
using special frequency converters is an important advantage 
of such a vibratory regime. A two-fold decrease of the vibra¬ 
tional frequency of the working unit with the amplitude re¬ 
maining unchanged decreases the transport velocity approxi¬ 
mately by the same factor. Consequently, to maintain the 
transport velocity the amplitude of the vibrations must be in¬ 
creased by a factor of 2, which does not present a problem 
for a resonance setup. The overall gain in the ratio of the 
strength and vibration insulation of the setup consists in a 
two-fold decrease of the force of inertia. In addition, the stiff¬ 
ness of the expensive elastic system decreases four-fold. 

StromizmeriteT CJSC is developing jointly with the 
Nizhny Novgorod State Technical University new high-ca¬ 
pacity energy-conserving vibratory conveyors and feeders 
for the glass and building materials industries on the basis of 
the low-frequency resonance vibratory mechanism proposed 
above. A resonance vibratory drive makes it easier to de¬ 
velop systems for automatically controlling technological 
processes and to attain high accuracy in dispensing and high 
stability in conveying loose and lumpy materials. 
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